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Abstract. Let A{n, d) be the maximum number of 0, 1 words of length n, any two having Hamming 
distance at least d. We prove ^1(20, 8) = 256, which implies that the quadruply shortened Golay code 
is optimal. Moreover, we show A(18,6) < 673, A(19,6) < 1237, A(20,6) < 2279, A(23,6) < 13674, 
A{19,8) < 135, A{25,8) < 5421, A(26,8) < 9275, A(21,10) < 47, A(22, 10) < 84, A(24, 10) < 268, 

O ■ ^(25, 10) < 466, A{26, 10) < 836, A{27, 10) < 1585, A{25, 12) < 55, and A{26, 12) < 96. 

OA ' The method is based on the positive semidefiniteness of matrices derived from quadruples of 

^.. . words. This can be put as constraint in a semidefinite program, whose optimum value is an upper 

C^ ' bound for A(n, d). The order of the matrices involved is huge. However, the semidefinite program is 

highly symmetric, by which its feasible region can be restricted to the algebra of matrices invariant 
under this symmetry. By block diagonalizing this algebra, the order of the matrices will be reduced 
so as to make the program solvable with semidefinite programming software in the above range of 
values of n and d. 
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For any k, we will identify elements {0, 1} with 0, 1 words of length k. The Hamming 
distance dH{v,w) between two words v,w is the number of i with Vi ^ Wi. 
Throughout we denote 

(1) iV:={0,ir. 



A code of length n is any subset C of iV. The minimum distance of a code C is the minimum 
*0 \ Hamming distance between any two distinct elements of C . Then A{n, d) denotes the 

^—N ' maximum size (= cardinality) of a code of length n with minimum distance at least d. 

Computing A{n, d) and finding upper and lower bounds for it have been long-time 
focuses in combinatorial coding theory (cf. Mac Williams and Sloane [11]). Classical is 
Delsarte's bound [3]. Its value can be described as the maximum A2{n,d) of 'Yliuv<^N ■^■u,v, 
where X is a symmetric, nonnegative, positive semidefinite N x N matrix with trace 1 
5^ I and with Xu,v = H u,v £ N are distinct and have distance less than d. Then A{n, d) < 

A2{n, d), since for any nonempty code C of minimum distance at least d, the matrix X with 
Xu V = \C\~^ ii u,v (z C and X^ ^ = otherwise, is a feasible solution with objective value 
\C\. 

This is the analytic definition of the Delsarte bound (in the vein of Lovasz [9], cf. [12], 
[16]). It is a semidefinite programming problem, but of huge dimensions (2"), which makes 
it hard to compute in this form. 
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However, the problem is highly symmetric. The group G of distance preserving permu- 
tations of the set A^ acts on the set of optimum solutions: if (Xu^v) is an optimum solution, 
then also (-'^7r(M),7r(t))) is an optimum solution for any it G G. Hence, as the set of optimum 
solutions is convex, by averaging we obtain a G-invariant optimum solution X. That is, 
^Tr{u),w{v) = Xu^y for all u, V and all n G G. So X^^v depends only on the Hamming distance 
of u and v, hence there are in fact at most n + 1 variables. Since (in this case) the algebra of 
G-invariant matrices is commutative, it implies that there is a unitary matrix U such that 
U*XU is a diagonal matrix for each G-invariant X. It reduces the semidefinite constraints 
of order 2" to 2" linear constraints, namely the nonnegativity of the diagonal elements. As 
the space of G-invariant matrices is n + 1-dimensional, there are in fact only n + 1 different 
linear constraints, hence it reduces to a small linear programming problem. 

So the Delsarte bound is initially a huge semidefinite program in variables associated 
with pairs and singletons of words in A^, that can be reduced to a small linear program, with 
a small number of variables. In [17] this method was generalized to semidefinite programs 
in variables associated with sets of words of size at most 3. In that case, the programs 
can be reduced by block diagonalization to a small semidefinite program, with a small 
number of variables. A reduction to a linear program does not work here, as in this case 
the corresponding algebra is not commutative. This however is not a real bottleneck, as 
like for linear programming there are efficient (polynomial-time) algorithms for semidefinite 
programming. 

In the present paper we extend this method to quadruples of words. Again, by a 
block diagonalization, the order of the size of the semidefinite programs is reduced from 
exponential size to polynomial size. We will give a more precise description of the method 
in Section [2j 

The reduced semidefinite programs still tend to get rather large, but yet for n up to 28 
and several values of d, we were able to solve the associated semidefinite programming up to 
(more than) enough precision, using the semidefinite programming algorithm SDPA. It gives 
the new upper bounds ^4(71, d) for A{n,d) displayed in Table 1. One exact value follows, 
namely ^4(20, 8) = 256. It means that the quadruply shortened Golay code is optimum. In 
the table we give also the values of the new bound where it does not improve the currently 
best known bound, as in many of such cases the new bound confirms or is very close to this 
best known bound. 

Since A{n, d) = A{n + 1, d + 1) if d is odd, we can restrict ourselves to d even. We refer 
to the websites maintained by Erik Agrell [1] and Andries Brouwer [2] for more background 
on the known upper and lower bounds displayed in the table. 

In the computations, the accuracy of the standard double precision version of SDPA 
(already considered in the comparison [13]) was insufficient for several of the cases solved 
here. The semidefinite programs generated appear to have rather thin feasible regions so 
that SDPA and the other high-quality but double precision codes terminate prematurely 
with large infeasibilities. We have used the multiple precision versions of SDPA developed 
by M. Nakata for quantum chemistry computations in [15]. Fortunately, the quadruple 
precision version was sufficient in all cases tabulated. The even higher precision versions as 
also implemented by the second author for interactive use at the NEOS Server [14] would 
have needed excessive computing times. Still, the times needed in Table 1 below varied 
from a few hours for the small cases to 1^/2 days for ^4(20, 4), 13 days for ^4(23, 6), 22 
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Table 1. Bounds for A{n,d) 



days for ^4(25,8), 30 days for ^4(27, 10), and 43 days for ^4(26,8). 

The approach outlined above of course suggests a hierarchy of upper bounds by con- 
sidering sets of words of size at most k, for k = 2,3,4, .. .. This connects to hierarchies 
of bounds for 0,1 programming problems developed by Lasserre [7], Laurent [8], Lovasz 
and Schrijver [10], and Sherali and Adams [19]. The novelty of the present paper lies in 
exploiting the symmetry and giving an explicit block diagonalization that will enable us to 
calculate the bounds. 

In fact, the relevance of the present paper might be three-fold. First, it may lie in coding 



and design theory, as we give new upper bounds for codes and show that the quadruply 
shortened Golay code is optimal. Second, the results may be of interest for algebraic com- 
binatorics (representations of the symmetric group and extensions), as we give an explicit 
block diagonalization of the centralizer algebra of groups acting on pairs of words from N. 
Third, the relevance may come from semidefinite programming theory and practice, by ex- 
ploiting symmetry and reducing sizes of programs, and by gaining insight into the border of 
what is possible with current-state semidefinite programming software, both as to problem 
size, precision, and computing time. 

We do not give explicitly all formulas in our description of the method, as they are some- 
times quite involved, rather it may serve as a manual to obtain an explicit implementation, 
which should be straightforward to derive. 

2. The bound Ak{n,d) 

For any n,d,k £ Z+, we define the number j4fc(n, d) as follows. Let M be the collection of 
codes 5" C A^ of minimum distance at least d. (Recall that N = {0, 1}".) For any t, let A/j 
be the collection of S* G AA with \S\ < t. 
For 5 G A/fc, define 

(2) M{S) := {S' e Af \ S <z S', \S\ + 2\S' \S\< k}. 

The rationale of this definition is that \S' U S"\ < k for all S", S" G Af{S). 

For X : A/fc -^ K and S G TVfc, let Ms{x) be the Af{S) x Af{S) matrix given by 



(3) Ms{x)s',s" 



for S',S" eMiS). Define 



x{S' U S") if S' U S" G AA, 
otherwise. 



(4) Ak{n,d) := max{N^ ^{{v}) \ x : A4 — )• M, x(0) = l,Ms{x) positive semidefinite 

veN 
for each 5 G A/^}. 

Note that, as x{S) occurs on the diagonal of Adsix), x has nonnegative values only. 
Proposition 1. A{n,d) < Ak{n,d). 

Proof. Let C be a maximum-size code of length n and minimum distance at least d. Define 
x{S) := 1 if 5 C C and x{S) := otherwise. Then Ms{x) is positive semidefinite for each 
S G J\fk. Moreover, A{n, d) = \C\ = E.^tv <M)- I 

The upper bound A2{n, d) can be proved to be equal to the Delsarte bound [3] (see [5]). 
The bound given in [17] is a slight sharpening of A^{n,d). 

Now to make the problem computationally tractable, let again G be the group of per- 
mutations of N that maintain distances. Then, if x is an optimum solution of (j4]) and 



TT € G, also x^ is an optimum solution. (We refer to Section [3121 for notation.) As the fea- 
sible region in ^ is convex, by averaging over all vr € G we obtain a G-invariant optimum 
solution. So we can reduce the feasible region to those x that are G-invariant. Then Mg(x) 
is Gg-invariant, where Gs is the G-stabilizer of S (= set of vr G G with tt{S) = S). This 
allows us to block diagonalize Ms{x), and to make the problems tractable for larger n. In 
the coming sections we will discuss how to obtain an explicit block diagonalization. 

It is of interest to remark that the equality A{20, 8) = 256 in fact follows if we take 
k = 4 and require in ^ only that Ms{x) is positive semidefinite for all S with \S\ = or 
\S\=A. 

An observation useful to note (but not used in the sequel) is the following. A well-known 
relation is A{n + 1, d) < 2A{n, d). The same relation holds for Ak{n, d): 

Proposition 2. For aUn,d: Ak{n + l,d) < 2Ajt{n,d). 

Proof. Let x attain the maximum ([!]) for Afc(n + l,d). For each S C {0, 1}", let S' := {wO \ 
weS} and S" := {wl \ w e S}. Define x'{S) := x{S') and x"{S) := x{S") for all S e Af. 
Then x' and x" are feasible solutions of ([!]) for Ak{n,d). Moreover J2veioi}"-(^'({'^}) + 
x"{{v})) = Y:^^{o,i}n+ix{{v}). Thns2Ak{n,d)>Ak{n + l,d). ' | 

This implies, using ^(20, 8) = 256 and yl(24, 8) = 4096 (the extended Golay code), that 
^4(21,8) = 512, ^4(22,8) = 1024, ^4(23, 8) = 2048, and ^4(24,8) = 4096. We did not 
display these values in the table, and we do not need to solve the corresponding semidefinite 
programming problems. 

From now on we will fix A; = 4, and we will use the name k for other purposes. In Section 
[6] we will discuss how to find a further reduction by considering words of even weights only, 
which is enough to obtain the bounds. 

3. Preliminaries 

In this section we fix some notation and recall a few basic facts. Underlying mathematical 
disciplines are representation theory and C*-algebra, but because the potential readership 
of this paper might possess diverse background, we give a brief elementary exposition. 

3.1. Notation 

We denote 

(5) [s,t] := {i £ Z \ s < i < t} and [s, t]cvcn := {i £ "^ \ s < i < t;i even}. 

Throughout this paper, P, T, and N denote the sets of ordered pairs, triples, and 
n-tuples of elements of {0, 1}, i.e., 

(6) P ■= {0, 1}2, T := {0, 1}3, and N := {0, 1}". 



As mentioned, we identify elements of {0, 1}* with 0, 1 words of length t. We will view {0, 1} 
as the field of two elements and add elements of P, T, and N modulo 2. 
For a : P — 7- Z+, we denote 

(7) ia :=a{W)+a{ll),ja := a(01) + Q(ll),nc, := a(00) + a(10) + a(01) +a(ll). 
For any finite set V and n E Z+, let 

(8) A^ :={A:y^Z+ I ^A(w) =n}. 

For any A G Ay, let 

(9) f^A :={/0: {1, •••,«} ^y I \p~'^iv)\ = A(v) for each t; G V}. 

So {il.x I A € Ay} is the collection of orbits on V^ under the natural action of the symmetric 
group Sn on F" (cf. Section El]) . 

Throughout, G denotes the group of distance preserving permutations of A^. The group 
consists of all permutations of coordinates followed by swapping and 1 in a subset of the 
coordinates. 

3.2. Group actions 

An action of a group H on a set Z is a group homomorphism from H into the group of 
permutations of Z. One then says that H acts on Z. An action of H on Z induces in a 
natural way actions of H on derived sets like Z x Z, V{Z), {0,1}^, and C^ . 

li TT £ H and z £ Z, then z^ denotes the image of z under the permutation associated 
with vr. If H acts on Z, an element z £ Z is called H -invariant if z'^ = z for each tt £ H. 
The set of iJ- invariant elements of Z is denoted by Z^ . 

A function (/> : Z — >• Z is H -equivariant if (piz^) = (f>{z)'^ for each z £ Z and each tt £ H. 
If Z is a vector space, the collection of i7-equivariant endomorphisms Z — )■ Z is denoted by 
Endiy(Z). It is called the centralizer algebra of the action of H on Z. 

If Z is a finite set and H acts on Z, then there is a natural isomorphism 

(10) EndH(C^) = (C^^^)^. 

If Z is a linear space, the symmetric group Sn acts naturally on the n-th tensor power Z*^". 
As usual, we denote the subspace of symmetric tensors by 

(11) Sym"(Z) := (Z®'^)^". 

3.3. Matrix *-algebras 

A matrix *-algebra is a set of matrices (all of the same order) that is a C-linear space and 
is closed under multiplication and under taking conjugate transpose (A i— )• X*). If a group 



H acts on a finite set Z, then [C^^^)^ is a matrix *-algebra. 

If A and B are matrix *-algebras, a function (/> : ^ — )• ;B is an algebra *-homomorphism 
if is linear and maintains multiplication and taking conjugate transpose. It is an algebra 
^-isomorphism if (p is moreover a bijection. 

li 4> : A ^ B is an algebra *-homomorphism and A £ A is positive semidefinite, then 
also 4>{A) is positive semidefinite. (This follows from the fact that any matrix X is positive 
semidefinite if and only if X = YY* for some linear combination Y of X, X'^ , . . . . This last 
statement can be proved by diagonalizing X.) 

The sets C™^'", for m G Z_|_, are the full matrix *-algebras. An algebra *-isomorphism 
A ^ B is called a full block diagonalization of ^ if ;B is a direct sum of full matrix *-algebras. 

Each matrix *-algebra has a full block diagonalization — we need them explicitly in 
order to perform the calculations for determining A4^{n,d). (A full block diagonalization is 
in fact unique, up to obvious transformations: reordering the terms in the sum, and resetting 
X I— 7- U*XU, for some fixed unitary matrix U, applied to some full matrix *-algebra.) 

3.4. Actions of ^2 

Let Z be a finite set on which the symmetric group 5*2 acts. This action induces an action 
of 5*2 on C . For it G {+, — }, let L± := {x G M | x'^ = itx}, where a is the non-identity 
element of 82- Then L+ and L„ are the eigenspaces of a. 

Let U± be a matrix whose columns form an orthonormal basis of L± . The matrices U± 
are easily obtained from the S'2-orbits on Z. Then the matrix [C/+ U-] is unitary. Moreover, 
U^XU^ = for each X in (C )^. As L+ and L_ are the eigenspaces of a, the function 
X ^ UXXU+ e UlXU^ defines a full block diagonahzation of {C^^^f^ 

3.5. Fully block diagonalizing Sym"(C^^^) 

We describe a full block diagonalization 

(12) X ■■ Sym"(C2x2) ^ clk,n-k]x[k,n-k]^ 

fe=0 

as can be derived from the work of Dunkl [4] (cf. Vallentin [20], Schrijver [17]). 
To this end, let, for any a G Ap, 

n 

(13) M„ := Y, (8)^P(0 e Sym"(C2x2). 

Here, for c = (ci,C2) G P, Ec denotes the {0,1} x {0,1} matrix with 1 in position ci,C2 
and elsewhere. Then {Ma \ a G Ap} is a basis of Sym"'(C^x^). (Throughout, we identify 
C2X2 with c^o^}x{o^},) 

For any a : P — > Z+ and k G Z+, define the following number: 



(14) 7a,fc := 

k 

E(_'i\V(k\ /m(01)+«(00)-fc\ /a(01)+a(ll)-fc\ /a(10)+a(00)-fc\ /«(10)+a(ll)-fc\\ 1/2 
V '-' \p)\\ a(01)-p A a(01)-p )\ a(10)-p )\ a(10)-p // ' 
p=0 

and the following [/c, n — /c] x [/c, n — /c] matrix r^^fc, where n := no,: 

, X ,„ X I 7q fe if i = ia and j = ja , 

(15) {T^kkj ■■={'''''' ^^ 

I U otherwise, 

for i,j G [/c, n — /c]. Now x is given by 

L|nJ 

(16) x(M„) = r„,fc 

k=0 

for a G A^. 

4. Fully block diagonalizing the matrices Ms{x) 

Recall that A/f consists of all codes of length n, of minimum distance at least d, and size 
at most t. Let x : Mi — )• M be G-invariant. For any S G ^4 and any tt & G, we have that 
Ms{x) is equal to M^(5)(x), up to renaming row and column indices. So we need to check 
positive semidefiniteness of Ms{x) for only one set S from any G-orbit on A/4. Moreover, 
Ms{x) belongs to the centralizer algebra of the G-stabilizer Gs of 5 (the set of all vr G G 
with 7r(S') = S). 

Consider any S G N4,. If \S\ = 4, then Ms{x) is a 1 x 1 matrix with entry x{S). So 
Ms{x) is positive semidefinite if and only if x{S) > 0. 

Moreover, if IS*! is odd and \S\ < 3, then Ms{x) is a principal submatrix of Mji{x), where 
R is any subset of S with \R\ = \S\ - 1. (This because if S' ^ S and |5| + 2\S' \ S\ < A, 
then \R\+2\S'\R\ < 4.) 

Concluding it remains to consider checking positive semidefiniteness of Ms{x) only for 
5 = and for one element from each G-orbit of codes SinM with l^j = 2. We first consider 
the case |5| = 2. 

Note that the G-orbit of any S G A/" with IS*] = 2 is determined by the distance m 
between the two elements of S. Hence we can assume S := {0, u}, where u is the element of 
N with precisely m I's, in positions 1, . . . ,m. Then there is a one-to-one relation between 

(17) iV' := {7; G iV I dniO, v) G [d, n] and dniu, v) G {0} U [d, n]} 

and M{S), given hy v >-^ S U {v}. 

For any v,w G N, let pi,^w '■ {!;••• 5^} — ^ P be defined by Pv,w{i) '■= {vi,Wi) for 
i = 1, . . . , n. This gives an embedding 



(18) $ : EndG's(C^^^)) ^ Sym'"(C2x2) ^ Sym"-™(C2x2) 

defined by 






(19) cD(X):= ^ ^5uW,5uwMp„,„ 

for XGEndGs(C-^(^)). 

The image of $ is equal to the hnear huh of those Mq, (g) M/3 with a e A^ and /? e Ap~™ 
such that ia + ip S [d, n], ja+j/3 £ [d, ra], m — ia + ip G {0}U [d, n], m — ja+jj3 G {0}U [d, n]. 

With the full block diagonalization (jl5|) this gives that the image is equal to the direct 
sum over A;, / of the linear hull of the submatrices of T^^k ® '^I3,i induced by the rows and 
columns indexed by («,?') with i + i' ^ [d, n] and m — i + i' € {0} U [d, n]. 

The stabilizer Gs contains a further symmetry, replacing any c G A^ by c + u (mod 2). 
This leaves S = {0, u} invariant. It means an action of 5*2, and the corresponding reduction 
can be obtained with the method of Section [3l 



5. Fully block diagonalizing ]\/[^{x) 

In this section we consider 5 = 0. Then M{S) = ^2, which is the set of all codes of length n, 
minimum distance at least d, and size at most 2. We are out for a full block diagonalization 
of the centralizer algebra EndG(C^) of the action of G on C^. This will be obtained in a 
number of steps. 

5.1. The algebra A. 

We first consider an algebra A consisting of (essentially) 4x4 matrices. For any c & P = 
{0, 1}^, let c := c + (1, 1) (mod 2). Let A be the centralizer algebra of the action of 5*2 on 
P generated by c 1— )■ c on c G P. We can find a full block diagonalization with the method 
of Section I3I4I We need it explicitly. Note that 

(20) A = {AeC^''^ \ A-^ = Ac,d for ah c,deP} 

and that ^ is a matrix *-algebra of dimension 8. 

For c,d ^ P, let E^.^^ be the 0, 1 matrix in C^^^ with precisely one 1, in position (c, d). 
Recall T = {0, l}^ and define for t G T: 

(21) Bt := E,^d + E--^, 

where (c, d) is any of the two pairs in P^ satisfying 

(22) ci + C2 = ti.di + d2= t2, C2 + d2 = is- 
Then {Bt | t G T} is a basis of A. 



For i £ {0, 1}, let U, G C^^{°'i> be defined by 

(23) mc,a = lV2{-ir^6a,c,+c, 

for c £ P and a G {0, 1}. Then one directly checks that the matrix U := [Uq Ui] is unitary, 
i.e., U*U = I. Moreover, for all c,d £ P and a, b, i,j £ {0, 1} we have 

(24) {U*E,,dUj)a,b = {Ui)c,a{Uj)d,b = \{-lt^^''''5a,c,+c,h,d.+d2- 

Hence, if i G T and c, d satisfy (j22p . then 

(25) {UtBtU.U = l^i-ir^^^"^ + (-l)^'=^+^''^^+^+^>a,e,+c,<5Ml+d2 = 

1 (_l)-2+.d.(i ^ (-l)*+:')5,,,,+,,5,,,,+,, = {-lf-'5ij5a,t,5i,,t,. 

So UqAUi = 0, and hence, as dim^ = 8, U*AU gives a full block diagonalization of A. 
Moreover 

(26) U*BtUi = {-lf'Et,,t2- 

5.2. The algebra Sym"(A) 

Our next step is to find a full block diagonalization of EndG(C ), where N"^ is (as usual) 
the collection of ordered pairs from A^ = {0, 1}". 
There is a natural algebra isomorphism 

(27) EndG(C^') ^ Sym"(^) 
by the natural isomorphisms 

(28) C(^°'^^")' ^ C^{0^}T ^ ((i-{0,l}2^®n^ 

using the fact that G consists of all permutations of A^ given by a permutation of the indices 
in {0, . . . , n} followed by swapping and 1 on a subset of it. 

Let Uq and Ui the P x {0, 1} matrices given in Section [5111 Define 



(29) 4> : Sym"(^) ^ Sym"^(C^°'i>^^°'^>) ® Sym"-'^(C^°'^>^^°'^>) 

by 

n 

(30) ct>{A) := (C/g®'" O [/^^"-'^)*^(^®'" ^ ^«5n-m^ 

m=0 

for A £ Sym"'(^). Trivially, (/> is hnear, and as U*AU = U^AUq © U^AUi, </> is a bijection 

10 



(cf. Lang [6], Chapter XVI, Proposition 8.2). Moreover, it is an algebra *-isomorphism, 
since U*Ui = I for i = 0,1 and hence U^"" O Uf""'"^ is unitary. 

Since a full block diagonalization of Syin™(C^^'^), expressed in the standard basis of 
Sym™(C^^^), is known for any m (Section I3l5p . and since the tensor product of full block 
diagonalizations is again a full block diagonalization, we readily obtain with (j) a full block 
diagonalization of Syni'*(^). To use it in computations, we need to describe it in terms 
of the standard basis of Sym"'(^). First we express (j) in terms of the standard bases of 
Sym"(^) and of Sym™(C2x2) and Sym"^™(C2x2). 

Let A^ and Qx be as in ([H]) and Q. For A G A^, define 

n 

(31) Ba:= 5](g)i?p«- 

Then {-Ba | A G A^;} is a basis of Sym''(^). 

We need the 'Krawtchouk polynomial': for n,k,t G Z_|_, 

k 

(32) KUt):=Y.i-^yQil-^- 

i=0 

For later purposes we note here that for all n, k, t: 

(33) Kll^.it) = (-l)*K«(t). 

For A G A^, a G A^, /3 G A^"'", define 

(34) T?A,a,/3 := 6x',a+p n ^a!c1)(/?(c)), 

cGP 

where for A G Ay, A' G Ap is defined by 

(35) A'(c) := A(cO) + A(cl) 

for c £ P. 

We now express (j) in the standard bases ([3T]) and ([13]). 

Proposition 3. For any A G A^, 

n 

(36) (/,(Sa) = Y1 ^x,a,pMa^M^. 

™=0a6A?J,/3eAf 



\.n — m 



Proof. By ()26p . the m-th component of (j){B\) is equal to 
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m n 

(37) E {®E^.(^),M^)) ® ( (2) i-^r^'^E^.(^lM^) 



lit II IIL 






If we sum over a := /i' and /3 := z^', we can next, for each c £ P, sum over j and set 
vicl) := J, and /i(cl) := A(cl) — j. In this way we get that the last expression in ([37|) is 
equal to 

A(cl) 

(38) E (n E(-i)^'(A(:!t.)('f))^°^^/' = 

QSA™,/3eAj-'" CgP j=0 
c«+/3=A' 



This describes the algebra isomorphism in (|29p . With the results given in Section I3I5I 
it implies a full block diagonalization 



(39) V':Sym"M)^0 d'^ 

described by 



m—k]x[k,m—k] ^ fri[l,n—ni—l]x[l,n—m—l] 
m=0 A:=0 1=0 



n LfH L|("-'Ti)J 

(40) V'(i?,) = Vn^,M(^A) 

m=0 fc=0 /=0 

where 

(41) ^m,kABx)-= E ^\,a,lS^a,k^Tfi,l 

aeA™,/3GAp,-'" 

for A S A^. Inserting (J14p and (J15p in (j4ip makes the block diagonalization explicit, and 
it can readily be programmed. Note that a, /3 in the summation can be restricted to those 
with a + (3 = A'. Note also that at most one entry of the matrix Ta^k ^ '^13,1 is nonzero. 

5.3. Deleting distances 

For m, k, I, we will use the natural isomorphism 
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(A'?\ fp{[k,m—k]x[l,n—m—l])x{[k,m-~k]x[l,n~m.~l]) r^ ^[k,m~k]x[k,m~k] ^ fp[l,n~m~l]x[l,n~m~l] 

Proposition 4. Let D C [0,n]. Then the linear hull of 

(43) ii^mMBx) I A G A^,ix',jx' G D} 
is equal to the subspace C^^^ of (j42p . where 

(44) F := {{i, i') (^ [k,m - k] x [l,n - m - l\\ i + i' £ D]. 



Proof. For any A G A^, if ijJm,k,l{B\)(^i^ii)^(^jj') is nonzero, then i + i' = i\i and j + j' = jy- 
This follows from ()40p and from the definition of the matrices F^ ^ (cf. (J15p ). 

Hence, for any fixed a,b £ Z4., the linear hull of the ipm,k,l{B\) with i\i = a and jy = b 
is equal to the the set of matrices in (|42p that are nonzero only in positions {i,i'), {j,j') 
with i + i' = a and j + j' = b. | 

So if distances are restricted to Z) C [0, n], we can reduce the block diagonalization to 
those rows and columns with index in F. 

5.4. Unordered pairs 

We now go over from ordered pairs to unordered pairs. First, let M2 '■= M2 \ {0}, and 
consider EndG(C 2). Let r be the permutation of iV^ swapping (c, d) and (d, c) in N"^ . Let 
Qr be the corresponding permutation matrix in C^ ^^ . Note that A^^ corresponds to the 
set of row indices of the matrices B\ (cf . ([25]) ) . 

For any A G A^, let A G A^ be given by A(ii, ^2, ^3) := A(ii, ^2, ^3 + ti) for t G T. So for 
any c G P, A(cl) = A(cl) if ci = and A(cl) = A'(c) - A(cl) if ci = 1. Then B^ = QrBx- 

Now (j33|) gives that for any m and a G Ap, /3 G Ap~™ one has 

(45) ^A,„,^ = (-l)'''^A,a,/3. 

This implies that the matrix iprn,k,l{B\ + Bj) has only O's in rows whose index («,«') has i' 
odd. Similarly, the matrix iprn,k,i{B\ — Br) has only O's in rows whose index {i,i') has i' 
even. As the space of matrices invariant under permuting the rows by r is spanned by the 
matrices Bx + B-^, this space corresponds under ipm,k,l to those matrices that have O's in 
rows whose index (i,i') has i' odd. 

A similar argument holds for columns if we consider A(ti,t2;i3) := -^(^1,^25^3 + ^2)- 
Hence the image of tpm,k,i of the set of elements of Sym"'(^) that are invariant under the 
operations A 1— )• A and A 1— )■ A is precisely equal to 



(461 fp[k,m-k]x[k,m-k] „ ip[l,n-m~l]cvcn'x[l.,n-m-l]c 
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5.5. Adding 

So far we have a full block decomposition of EndcCC-'^), We need to incorporate in it. It 
is a basic fact from representation theory that if Vi, . . . , Vf is the canonical decomposition 
of C-'^ into isotypic components (cf. Serre [18]), then EndG(C^2) = 0j=i EndG'(yj), and 
each EndG(Vi) is *-isomorphic to a full matrix algebra. 

We can assume that Vi is the set of ff-invariant elements of C'^. Hence, as is G- 
invariant, V( := C® © Vi is the set of G- invariant elements of C-^. One may check that 
the block indexed by {m,k,l) = (n, 0, 0) corresponds to Vi. So replacing block (n,0,0) by 
EndG(^i') gives a full block diagonahzation of EndclC-^"). Note that EndciVl) = End(y/). 

We can easily determine a basis for V(, namely the set of characteristic vectors of the 
G-orbits of A2. Then for any B G EndG'(C-'^), we can directly calculate its projection in 
Eind{V(). This gives the required new component of the full block diagonahzation. 

6. Restriction to even words 

We can obtain a further reduction by restriction to the collection E of words in {0, 1}" of 
even weight. (The weight of a word is the number of I's in it.) By a parity check argument 
one knows that for even d the bound A{n, d) is attained by a code C Q E. A similar 
phenomenon applies to Ak{n,d): 

Proposition 5. For even d > 2, the maximum value in ^ does not change if x{S) is 
required to be zero if S % E. 

Proof. Let e : N ^ E he, defined by £{w) = w \i w has even weight and £{w) = it; + e^ if 
w has odd weight. Here e„ is the n-th unit basis vector, and addition is modulo 2. If d is 
even, then for all v,w £ N: dH{v,w) > d \i and only if dH{£{y),e{w)) > d. Now e induces 
a projection p : M.-^ — t- M , where £ is the collection of codes in M with all words having 
even weight. 

One easily checks that if Ms{x) is positive semidefinite for all S, then Ms{p{x)) is 
positive semidefinite for all S. Moreover, Ylv&N Pi^)i{'^}) ~ S^jsat ^({^D- I 

This implies that restricting x to be nonzero only on subsets S of E does not change 
the value of the upper bound. However, it gives a computational reduction. This can 
be obtained by using Proposition [4] and by observing that the restriction amounts to an 
invariance under an action of S2, for which we can use Section [3141 The latter essentially 
implies that in (j46p we can restrict the left hand side factor to rows and columns with index 
in [k,m — A:]even- As it means a reduction of the program size by only a linear factor, we 
leave the details to the reader. 

Acknowledgement. We thank Niels Oosterling for very helpful comments on the method. 
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